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Abstract 

Using an isospin-dependent quantum molecular dynamics (IQMD) model, nuclear stopping is 
analyzed in asymmetric colliding channels by keeping the total mass fixed. The calculations have 
been carried by varying the asymmetry of the colliding pairs with different neutron-proton ratios in 
center of mass energy 250 MeV/nucleon and by switching off the effect of Coulomb interactions. We 
find sizable effect of asymmetry of colliding pairs on the stopping and therefore on the equilibrium 
reached in a reaction. 
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I. INTRODUCTION 



In recent years, the study of heavy-ion colhsions from low to relativistic energies has 
been focused on variety of phenomena that includes the multi- fragmentation [l| , anisotropy 
in the momentum distribution , as well as global stopping of the nuclear matter 
As we know, nuclear stopping is one of the essential observables that depends crucially on 
the reaction dynamics. The stopping has also been linked with the degree of thermalization 
and equilibrium reached in a reaction. This problem has been handled both theoretically 
and experimentally in recent years. Some have also tried to correlate it with the production 
of light charged particles. Even the role of symmetry energy has also been explored in this 
domain. It still remains to be seen how stopping (and/or thermalization/equilibrium) is 
affected when the asymmetry of the reacting partners is altered. As noted, the dynamics 
and energy deposition in asymmetric reaction can be quite different than in a symmetric 
reaction. The asymmetry of a reaction has also been reported to affect the collective flow, 
balance energy as well as elliptical flow of the colliding pairs. It is worth mentioning that 
the outcome and physical mechanism behind the symmetric and asymmetric reactions are 
entirely different. 

Following the establishment of radioactive beam facilities in many laboratories, it became 
possible to study the neutron-rich (or proton-rich) nuclear collisions at intermediate 
energies. The idea of studying nuclear stopping phenomena was introduced by Bass 
et al, jsl via the 'isospin-mixing' method. In 1988, Bauer 0], pointed out that the 
nuclear stopping at intermediate energies is determined by the mean-field as well as by 
the in-medium NN cross sections. In 1998, Li et at, jol] found that the degree of isospin 
equilibrium depends sensitively on both the in-medium NN cross section and equation of 
state of asymmetric nuclear matter. Another study in 2001 [3], explored the possibility of 
using nuclear stopping to probe the isospin dependence of in-medium NN cross-section. In 
2002 [sl, authors studied the behavior of excitation function Q^z/nudeon and concluded 
that Qzz/nucleon can provide information about the isospin dependence in terms of 
binary cross-sections. The recent work of many authors suggested that the degree 

of approaching isospin equilibration helps in probing the nuclear stopping in HIC. Several 
more studies also focused in recent years on the isospin degree of freedom In 2006, one 
of us and co-workers [lo| correlated the multifragmentation with global nuclear stopping. 
Their study revealed that the light charged particles (LCP's) act in a similar manner as the 
anisotropy ratio. They, however, excluded the isospin content of the colliding nuclei. In a 
recent communication jll|, one of us and co-workers tried to study the effect of symmetry 
energy and isospin-dependent cross-section on nuclear stopping. Our findings revealed 
that the degree of stopping depends weakly on the symmetry energy and strongly on the 
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isospin-dependent cross-section. Most of the above mentioned theoretical and experimental 
studies concentrated on the dynamics for symmetric reacting partners. 
It, however, remains to be seen how asymmetry of a colliding pair affects the stopping 
and thermalization. As noted by FOPI collaboration [l^, asymmetry of colliding pairs 



13| can play decisive role in reaction dynamics. We plan to address this question in this 



present paper. In order to focus exclusively on the asymmetric aspects, we shall vary 
the masses of the projectile and target in such a manner that total system mass remains 
constant. Another hidden parameter that can affect outcome in asymmetric colliding nuclei 
is incident energy, we shall also keep center of mass energy fixed in present study. This 
study is conducted within the framework of IQMD model which is explained in the Sec.-II. 
The results are presented in Sec. -Ill, leading to the conclusions in Sec.-IV. 



II. THE MODEL 

The isospin-dependent quantum molecular dynamics (IQMD)[l^ model treats different 
charge states of nucleons, deltas and pions explicitly, as inherited from the Vlasov-Uehling- 
Uhlenbeck (VUU) model. The IQMD model was used successfully in analyzing the large 
number of observables from low to relativistic energies. The isospin degree of freedom 
enters into the calculations via both cross sections and mean field. 
In this model, baryons are represented by Gaussian-shaped density distributions 

„, , 1 -(r-r,Xt))2 -(p-p, (t)9 .2L 

fi{r,p,t) = ——e 2^ e ^ . (1) 

Nucleons are initialized in a sphere with radius R = 1.12 A^^^ fm, in accordance with the 
liquid drop model. The Gaussian width which we are taking is in close agreement as 
given in ref. jl5|. Each nucleon occupies a volume of so that phase space is uniformly 
filled. The initial momenta are randomly chosen between and Fermi momentum pp. 
The nucleons of the target and projectile interact via two and three-body Skyrme forces 
and Yukawa potential. The isospin degrees of freedom is treated explicitly by employing 
a symmetry potential and explicit Coulomb forces between protons of the colliding target 
and projectile. This helps in achieving the correct distribution of protons and neutrons 
within the nucleus. 

The hadrons propagate using Hamilton equations of motion: 

dfi d < H > dpi d < H > 

dt dpi ' dt dfi 
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with 

< H >=< T > + < > is the Hamiltonian. 

^2 



+ mp,t)V''ir\r)f^ir',p',t)dfdr'dpdp'. (3) 



The baryon-baryon potential V^^ , in the above relation, reads as 



^ y ') — ^ Skyrme ' ^Yukawa ' ^Coul ' ^ Sym 

= ti6{r' - f) + t26{r' - f)p^-\^-^] 

exp{\ r' -f\/fi) ZiZje^ 
+ ^3 — — r 



(I r' — r \r' — r\ 

+ u]-TlTlb{r[~r-). (4) 

Po 

Where /i = 1.5/m, ^3 = — 6.66MeV^, ^4 = lOOMeV^. The values of t\ and ^2 depends on the 



values of a, /3, and 7 [3[. Here Zj and denote the charges of the i^^ and j*^ baryon, and 
T3, T| are their respective T3 components (i.e. 1/2 for protons and -1/2 for neutrons). The 
parameters \i and ti, , are adjusted to the real part of the nucleonic optical potential. 



III. RESULTS AND DISCUSSION 

In the present study, projectile mass is varied between 16 and 56 units and targets 
are chosen as different isotopes of Xe, S'n, Ru in such a way that total mass of the 
reaction remains constant (= 152) for all channels. For example, we take the reactions of 
80^6 xei36, 14^^28 +54 Xei24, .g^^^ +50 Sn^''\ 2oCa^' +50 Sn'^^ ^aCt^^ +44 Ru^''\ and 
2%Fe>^ +44 Ru^^ etc. Although, the total mass remains constant, the asymmetry of the 
reaction 77 = | {At — Ap)/{At + Ap) \ keeps varying between 0.2 and 0.7. Nuclear stopping 
in HIC has been studied with the help of different variables. A direct measure of nuclear 
stopping is the rapidity distribution defined as [3] 

^W-2^>,)-pil(.)' 

where E{i) and Pz{i) are respectively, the energy and longitudinal momentum of the i^^ 
particle. For a complete stopping, one expects a single peaked Gaussian. Obviously, narrow 
Gaussian indicates better thermalization (equilibrium) compared to broader ones. 
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FIG. 1: The rapidity distribution dN/dY as a function of reduced rapidity for free nucleons, LMF's 
and IMF's respectively at the center of mass energy Eq.m. = 250 MeV/nucleon. Different curves 
correspond to different asymmetries varying from 0.2 to 0.7 for semi-central impact parameter. 



Another quantity used in the hterature [l6| is anisotropy ratio (R) defined as 



U^'-§^ (6) 

where summation runs over all nucleons. The transverse p±{i) and longitudnal p\\{i) mo- 
menta reads respectively as ■\Jpl{i) + Py{i) and Pz{i)- Naturally, for a complete stopping R 
should be close to unity. Some studies use quadrupole moment Q^z to analyze the stopping 
and thermalization. Quadrupole moment Qzz is defined as 

Qzz = Y.^2p\{^)-pl{^)] (7) 

i 

Naturally, for a complete stopping, Q^z should be close to zero. 

Fig.l shows the rapidity distribution ^ for the emission of free nucleons as well 
as light mass fragments (LMF's) [(2 < A < 4)], and intermediate mass fragments 
(IMF's) [(5 < A< Atot/Q)] at a fixed centre-of-mass energy. Note that the asym- 
metry varies between 0.2 and 0.7 corresponding to the reactions of sO^^ +54 Xe^^^ 
{t] = 0.7), uSz^' +54 Xei24 = 0.6), ig^^' +50 ^n^^o ^ g.S), ^oCa^' +50 Sn'^' {r, = 0.4), 
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24Cr^° +44 Rv}^"^ (77 = 0.3), 2(yFe>^ +44 Ru^^ (77 = 0.2). The Yc.M./Yi,eam = corresponds to 
mid-rapidity (participant) zone and hence is responsible for the hot and compressed zone. 
On the other hand, Yc.M./Ybeam 7^ corresponds to spectator zone iXc.M./Yheam < ~1 
corresponds to target hke (TL) and Yc.M./Ybeam > 1 corresponds to projectile hke (PL) 
distributions). Interestingly, reaction corresponding to 77 = 0.7 yields peak at mid-rapidity 
which shifts towards negative side when one considers nearly symmetric reactions. We see 
that majority of free particles and LMF's are emitted from the mid-rapidity region, whereas 
IMF's are emitted from the spectator matter. Prom the shape of the Gaussian, one sees 
that free particles and LMP's are better indicator for the thermal source at 77 = 0.7. 
If these reaction channels are analyzed at a fixed lab energy, the situation would have 
been entirely different. In that case, nearly symmetric reactions 77 = 0.2 would yield better 
thermalized source. It is worth mentioning that in many studies, one has kept the lab 
energy fixed. The absolute conclusion in these studies are not due to asymmetry alone. 

In Pig. 2, we display the impact parameter dependence of global variables (R and 
Qzz) along with the multiplicity of LMP's. The results are displayed at Eq.m. — 250 
MeV/nucleon. Different curves in each panal represent different asymmetries. We observe 
that R and 1/Qzz behave in a similar fashion (note that R and Q^z will behave in just 
opposite fashion). The amount of stopping/equilibrium increases with the asymmetry of 
the reaction. As we know, major contribution for the stopping of nuclear matter comes 
from the hot and compressed matter that decreases almost linearly with impact parameter. 
To correlate the degree of stopping with the multiphcity of light charged particles, we show 
in the last panal the impact parameter dependence of the multiplicity of light charged 
particles. The behavior of light charged particle with impact parameter is similar to that 
of anisotropic ratio and inverse of quadrupole moment. The decrease in the multiplicity of 
LMP's complements with corresponding increase in the heavy fragments. These fragments 
are the remanant of the spectator matter. Therefore, a decrease in the multiplicity of LMF's 
with impact parameter measures directly the decrease in the degree of equilibrium and 
hence global stopping. This implies that the LMP's production can act as an indicator for 
nuclear stopping. At the same time, we also note that behavior with respect to asymmetry 
does behave in same fashion. This happens due to the fact that participant zone increases 
when one moves towards nearly symmetric reactions. Therefore, more LMP's are produced 
for Tj = 0.2 compared to rj = 0.7. It is worth mentioning that in the case of multiplicity 
of LMP's, we talk about the number of particles while in the case of stopping, it is the 
momentum phase space (x, y and z components of momentum) that is responsible for the 
growth. 
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FIG. 2: The anisotropy ratio R, quadrupole moment 1/Qzz and multiplicity of LMF's as a function 
of impact parameters. Different curves correspond to different asymmetries varying from 0.2 to 0.7 
for semi-central impact parameters. 

The stopping at any time during the collision can be divided into the contributions 
emerging from the protons and neutrons. We decompose the stopping into the contributions 
due to protons and neutrons. Here, at each time step during the collision, stopping due to 
neutrons and protons is analyzed separately. 

Fig. 3(a) shows the final state quadrupole moment l/Qzz decomposed into contributions 
due to neutrons and protons as a function of asymmetry 77. Fig.3(b) shows the final state 
anisotropy ratio < i? > as a function of the asymmetry of the system. We see no difference 
between the contributions due to neutrons and protons. This is due to the fact that < R> 
is the ratio of the mean transverse momentum p±{i) to the mean longitudnal momentum 
p\\{i) —Pzii)- To see the clear contribution of neutrons and protons, one has to look into 
the contributions of transverse and longitudnal momenta as shown in Fig. 3(c) and 3(d), 
respectively. A linear enhancement in the transverse and longitudnal momentum can be 
seen with asymmetry 77. Further the contribution of neutrons exceeds the corresponding 
contributions due protons. 

It will be of further interest to see whether the above findings depend on the isospin 
asymmetry (N/Z dependence) or not. For this, we display in Fig.4, the anisotropy ratio 
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FIG. 3: The quadrupole moment \/Qzz and anisotropy ratio R as a function of asymmetry r] 
at a center of mass energy Ec.m. = 250 MeV/nucleon. The third and fourth panals show the 
variation of transverse and longitudnal momentum with asymmetry, respectively. Different curves 
correspond to the contribution of neutrons and protons along with total contribution. 

< R>, inverse of quadrupole moment Qzz and multiplicity of LMF's as a function of neutron 
to proton ratio (N/Z) at different impact parameters ranging from central to peripheral one. 
We are also showing the results at impact parameter 6 = 0.1 by changing the Gaussian 



width according to the formulae a = 0.16A^^/'^ + 0.49 as given in ref ISj]. Interestingly, 
we see no change in the results by the variation of Gaussian width. An increase in the 
number of neutrons will increase the number of collisions and hence the absolute value 
of < i? >, 1/Qzz will increase with N/Z ratio. Our exclusive findings are: (1) maximum 
stopping is obtained for the systems with larger neutron content. This is true at all colliding 
geometries. This dependence diminishes as one moves from central to peripheral geometry. 
It is due to the fact that nuclear stopping is governed by the participant zone only. Moreover, 
in systems with more neutron content, role of symmetry energy could be larger, whereas 
effects due to isospin-dependent cross-section could play a dominant role in systems with less 
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FIG. 4: (Color online) Same as Fig. 2, but as a function of N/Z ratio. Different curves correspond 
to different impact parameters ranging between central to peripheral one. 

neutron content (proton content). Therefore there is a possibility that the combined effect 
of symmetry energy and cross-section could be approximately the same for all systems with 
different neutron and proton content |12]. Also neutron- neutron or proton-proton cross- 
section is a factor of 3 lower than the neutron-proton cross-section. (2) On the other hand, 
the multiplicity of LMF's follows the reverse trend. This is due to the reason that in case of 
N/Z = 1.1 [rj = 0.2), participant zone is more compared to N/Z = 1.4 [r] = 0.7). Therefore, 
more LMF's are produced for N/Z = 1.1 as compared to N/Z = 1.4. 



IV. CONCLUSION 

Using the isospin-dependent quantum molecular dynamics (IQMD), we have studied the 
nuclear stopping in asymmetric colliding channels by keeping total mass fixed. The calcu- 
lations have been carried out by varying the asymmetry of the colliding pairs with different 
neutron-proton ratios in the center of mass energy 250 MeV/nucleon and by switching off 
the effect of Coulomb interactions. The contribution of the neutrons and protons is checked 
in terms of anisotropy ratio < R> and quadrupole moment Qzz- The maximum stopping 
is obtained for the systems having maximum neutron-proton ratio because the contribution 
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of neutrons remains enhanced throughout the asymmetry range. Moreover, this dependence 
becomes weaker as as one moves from central to peripheral geometry. Interestingly, reverse 
trend is obtained when we vary the multiphcity of LMF's with N/Z ratio. 
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